Border tissue morphology is associated with the pattern of visual field progression in open-angle glaucoma

The etiology of open-angle glaucoma (OAG) is yet unclear. This study investigated possible risk factors, such as the morphology of the border tissue that affect the pattern of visual field (VF) progression in eyes with OAG. 166 eyes of 166 OAG patients with an externally oblique border tissue (EOBT) at least in one direction were included. EOBT was obtained by analyzing enhanced depth imaging spectral-domain optical coherence tomography images. A pointwise linear regression was used to determine VF progression by measuring the deterioration rate of each point in the VF. The odds ratio of VF progression for each risk factor was estimated using logistic regression analysis. Seventy (42.2%) eyes showed VF deterioration. In multivariate analysis, longer follow-up period, higher baseline intraocular pressure (IOP), lower mean ocular perfusion pressure (MOPP), and smaller angular location of the longest EOBT were associated with VF progression (all p values were below 0.05). In the multivariate analysis, the location of the longest EOBT was significantly associated with inferior (p = 0.002) and central (p = 0.017) VF progression. In conclusion, VF progression pattern in OAG eyes is associated with the location of the longest EOBT as well as other known risk factors.

Factors associated with visual field progression locations. In univariate analysis, only a longer follow-up period was associated with the presence of overall VF progression in OAG eyes (p < 0.001). However, Table 1. Participant demographics and clinical characteristics. DBP diastolic blood pressure, EOBT externally oblique border tissue, IOP intraocular pressure, MD mean deviation, MOPP mean ocular perfusion pressure, PSD pattern standard deviation, SBP systolic blood pressure, VF visual field. Data are presented as the means ± standard deviations unless otherwise indicated.
A risk factor associated with VF progression in the inferior hemifield in univariate analyses was longer followup period (p = 0.0002). In multivariate analysis, the risk factors were longer follow-up period, worse baseline MD, and smaller longest EOBT angular location (p = 0.0004; p = 0.039; p = 0.005, respectively) ( Table 5).
Risk factors associated with central VF progression were longer follow-up period and greater baseline IOP in univariate analysis (p = 0.001; p = 0.038, respectively). In multivariate analysis, central VF progression was associated with a longer follow-up period, lower mean MOPP, and smaller longest EOBT angular location (p = 0.0001; p = 0.01; p = 0.009, respectively) ( Table 6).

Discussion
The present study demonstrated that the location of VF progression was significantly associated with the longest EOBT location in OAG eyes (Fig. 2). This result confirms the previous finding that EOBT has a topographical association with glaucomatous damage 14,15 , and further identifies preferential sites of VF progression associated with detailed deep ONH morphology. In addition, other factors including mean MOPP, baseline IOP, and mean IOP also affected specific VF progression patterns in OAG eyes. The present study is the first to demonstrate a possible relationship between the longest EOBT and VF progression patterns in OAG eyes over the long-term period. Figure 1 clearly shows the main finding of the present study, namely that VF progression pattern varies according to the location of the longest EOBT. Moreover, this relationship was more pronounced when the location of the longest EOBT was far from the macula. If the longest EOBT location was over + 45 degrees (group 5), VF progression occurred mostly in the superior hemifield. In contrast, if the longest EOBT location was at less than -45 degrees (group 1), the initial VF defect was inferiorly dominant at baseline, and further VF progression occurred more in the inferior hemifield than other groups. In other groups, the initial VF defect pattern was relatively superiorly dominant, but further overall VF progression appeared to occur simultaneously in superior and inferior VFs.
In previous studies, relationships between ONH morphology and preferential location of VF progression have been detected 14,[17][18][19][20][21] . These previous studies used parameters such as parapapillary atrophy (PPA) or disc tilt as biomarkers indicative of ONH morphology. In the present study, EOBT measured by OCT was used as a biomarker of deformation around the ONH 22 . EOBTs were found frequently, even in ONH without PPA; www.nature.com/scientificreports/  www.nature.com/scientificreports/ consistently, the longest EOBT was significantly associated with the VF defect pattern in non-myopic eyes as well as myopic eyes in a previous study 14 .
Nevertheless, the structural characteristics of the EOBT could not completely explain VF progression patterns in OAG patients. For example, VF progression in the superior field was not significantly associated with longest EOBT location. There may be several possible reasons for this. First, most patients already had a damaged superior hemifield defect initially except those in group 1, and further VF progression could have been less substantial in follow-up exams. However, the initial VF defect pattern of each group was associated with ONH deformation, consistent with previous studies 14,15,21 . Second, factors other than EOBT could contribute simultaneously to the VF progression pattern. In the present study, a thinner CCT and lower mean MOPP were associated with superior VF progression. Thinner CCT [23][24][25][26] and lower mean MOPP 9,27-29 are well known progression risk factors in glaucoma. Worse baseline MD was only associated with inferior VF progression, but not with superior VF progression in the OAG eyes. Considering that superior hemifield was mostly damaged to some extent initially especially when VF had worse MD, a greater chance of VF progression might occur relatively in the inferior VF in the eyes with longer axial length and worse baseline MD.
Lower mean MOPP was associated with central VF progression as well as superior VF progression in the current study. This result is consistent with previous reports that vascular factors including hypotension,  www.nature.com/scientificreports/ migraine, Raynaud phenomenon, blood pressure, and unstable MOPP are associated with central VF defects or progression 8,9,27 . However, no previous study has reported a relationship between lower MOPP and VF progression in the superior hemifield. In previous studies, changes in vascular structures located in superficial and deep ONH layers were usually found in relatively late-stage glaucoma 16,30,31 . Considering that superior VFs consistent with an inferior RNFL defect were present at initial diagnosis in the current study, corresponding vascular structures such as superficial and deep capillaries around the ONH in that area could be damaged more than other areas. We speculate that this could induce more significant VF aggravation if there is lower blood perfusion to those damaged capillaries. Central VF progression was also associated with higher baseline IOP in OAG eyes. In the previous study, higher IOP was associated with VF progression in overall VF including central 10° hemifield in normal-tension glaucoma 32 .
Meanwhile, we speculated that lower mean IOP might be the result instead of the cause of the overall VF progression. If the eyes showed VF progression, the clinician might add more IOP-lowering medications and eventually induce lower mean IOPs. For the limitation of the study design, this can be hardly proven for now. Further studies to show the causal relations between the factors are warranted.
Glaucoma is defined as a progressive optic neuropathy with unique VF patterns such as nasal step or arcuate scotoma that are usually located in the superior hemifield. The VF pattern is associated with glaucomatous RNFL defects of ONH. Previous studies have shown that geometric features around the ONH such as disc tilt or PPA are associated with initial VF defect patterns in glaucoma eyes with myopia 11,12 . Considering that myopic changes associated with ONH may be related to axial elongation during the growth period, regional structural differences around the optic nerve may also exist in eyes with various axial lengths 22 . Different ONH regional features such as EOBT might confer different susceptibilities of structures around the ONH to the given IOP in the eye. We hypothesized that IOP-related stress and strain around the ONH may induce initial structural damage including lamina cribrosa and RNFL, especially at the most vulnerable region, which corresponds to the longest EOBT. In OAG eyes, different EOBTs are associated with structurally weak areas around the ONH and induce initial glaucomatous damage and further VF progression for a given IOP in tandem with known risk factors such as ocular perfusion pressure.
There are several limitations to this study. First, we evaluated a relatively small number of cases. Thus, in several sub-analyses, some factors were close to but did not reach statistical significance. For compensating this limitation, we analyzed the factors associated with VF progression "including" specific patterns and "only with" specific patterns and tried to find the possible associated factors as many as possible. Second, the present study was designed as a retrospective longitudinal study. Because of the limitations of this study design, it was not possible to determine causal relationships between factors. For example, associated factors such as lower mean IOP can be the cause or the result of the VF progression in central VF. Third, OAG eyes at different stages were included in the analysis. If some region in the VF was already involved to some extent, there may have been less progression in that region. Fourth, The OAG eyes without EOBT were not included in this study. Although, this may arise selection bias, the reason for selecting such group was because an eye without EOBT results in a missing data. The main purpose of this study was to see the correlation between the location of the maximum EOBT and direction of progression of visual field defect. Fifth, the definition of visual field progression used in this study is different from the conventional definition. 42.2% of the patients showed VF progression in this study, which is significantly higher than the literature. If the existing definition is a macroscopic point of view, the definition used in this study was a microscopic point of view that sees changes from the point of view of each point in the Table 6. Risk factors associated with visual field progression in the central visual field. CI confidence interval, EOBT externally oblique border tissue, IOP intraocular pressure, MD mean deviation, MOPP mean ocular perfusion pressure. The best multiple logistic regression model was chosen based on the Akaike Information Criterion. Factors with statistical significance are marked in bold. www.nature.com/scientificreports/ visual field map. Therefore, visual field examination findings that were not clinically determined to be progression were considered progression in this study. Sixth, there is an overlap between the subgroups in subgroup www.nature.com/scientificreports/ analyses such as the central progression group and the superior or the inferior progression groups. However, the reason for this classification is that, the main objective of this study was to compare between the group with VF progression only in the superior VF and the group with VF progression only in the inferior VF. The secondary objective was to find the risk factors of the central VF progression. The central VF defect was further analyzed because it is crucial in clinical practice in view of the fact that it can affect the visual acuity.

Risk factors
In conclusion, the location of the longest EOBT and other factors such as mean MOPP and mean IOP contribute to specific VF progression patterns in OAG eyes.

Methods
This study was retrospective in nature, and subjects comprised OAG patients who visited Samsung Medical Center (Seoul, Korea) for at least 5-years of follow-up between July 2007 and December 2020. This study followed all guidelines for experimental investigation in human subjects, was approved by the Samsung Medical Center Institutional Review Board (2020-07-141-001), and adhered to the tenets of the Declaration of Helsinki. Informed consent was waived by the Samsung Medical Center Institutional Review Board.
Inclusion criteria were a diagnosis of OAG in addition to (1) the existence of externally oblique border tissue on enhanced depth imaging (EDI) spectral domain (SD)-OCT scans and (2) mean deviation ≥ − 12 dB on the first VF test. Exclusion criteria were concomitant ocular or systemic diseases that could affect VF tests, such as a history of vision-threatening retinal disease (e.g., retinal detachment, retinal vein occlusion) or neurologic disease. If glaucoma existed bilaterally, one of the two eyes was randomly included in the study.
Each participant underwent a comprehensive ophthalmic examination, including slit-lamp biomicroscopy, Goldmann applanation tonometry, gonioscopic examination, dilated stereoscopic examination of the ONH, color and red-free fundus photography (TRC-50DX model; Topcon Medical System, Inc., Oakland, NJ, USA), automated perimetry using a central 30-2 Humphrey field analyzer (HFA model 740; Humphrey Instruments, Inc., San Leandro, CA, USA) with the Swedish interactive threshold algorithm standard, axial length measurement (IOL Master®; Carl Zeiss Meditec, Jena, Germany), and ultrasonographic pachymetry (Tomey SP-3000; Tomey Ltd., Nagoya, Japan). The extent of the VF defect was measured using MD and PSD. VF tests were repeated once or twice a year, depending on the patient's condition. Reliable VF analysis was defined as a false-negative rate of < 15%, a false-positive rate of < 15%, and a fixation loss of < 20%. Unreliable VF test results were excluded from the analyses.
IOP at baseline was defined as the IOP at the first visit. The number of IOP-lowering medications at first visit was counted according to type (e.g., if a patient used a fixed combination, the patient was considered to be taking two medications). SBP and DBP were measured once at the first visit. MOPP was obtained using the formula MOPP = 2/3 [DBP + 1/3 (SBP-DBP)]-IOP.
Classification of progression pattern of the visual field defect. For each patient, we calculated the rate of VF loss (decibels per year) at each test location of the total deviation plot. Test locations were categorized into 10 zones according to the normal anatomy of the retinal nerve fiber layer 33 . Test points that can show high variability due to lens rim artifacts or proximity to physiological scotoma were excluded. Therefore, a total of 44 points were included in the study (Fig. 3). We considered a rate of deterioration faster than − 1 dB/year with sensitivity at the 1% level in any one point in zones 1 and 6 or any two points in the other eight zones as clinically significant deterioration based on the modification of the result of the previous report 34 . If significant deterioration was detected in zones 1 and 6, it was classified as central VF progression. Superior VF progression was defined when the VF progression pattern included significant deterioration in zones 1 to 5, and inferior VF progression was defined when the VF progression pattern included significant deterioration in zones 6 to 10. For clear comparison, 16 cases with progression in both superior and inferior zones were excluded. Measurement of the length and location of the longest externally oblique border tissue. Measurement of EOBT length by EDI SD-OCT was performed on each patient according to methods described in detail in a previous study 16 . Forty-eight radial-line B-scans (each at an angle of 3.75°) centered on the optic disc were obtained and every other section (each at an angle of 7.5°) was used. We defined the length of the EOBT as the distance between the inner end point of the EOBT and Bruch's membrane opening (Fig. 3). If the OCT image quality was too poor to recognize BMO or border tissue, the next scan section was used. The length of EOBT was measured in each scan and the scan with the longest EOBT was selected. A fundus photo of the same eye was superimposed on the OCT scan to identify the location of the fovea. The line connecting the center of the BMO and the fovea was defined as the fovea-BMO axis. The angle between the radial line and the fovea-BMO axis was defined as the angular location of the maximum EOBT. If a radial line was located below the fovea-BMO axis, it was regarded as a positive angle. Otherwise, it was regarded as negative angle.
Two independent observers (HJK and EDS) measured the length and location of the longest EOBT independently. The average length and location of the longest EOBT were used in subsequent analyses.
Statistical analysis. Statistical analyses were performed using R software v3. 6 www.nature.com/scientificreports/ were against the group without VF progression. In multivariate analysis, stepwise variable selection procedures based on the Akaike Information Criterion were implemented. ORs are presented as means with 95% confidence intervals. Two-sided P values are reported for all tests; P values less than 0.05 were considered statistically significant.  www.nature.com/scientificreports/

Data availability
The datasets that were analyzed for this study are available from the corresponding author on a reasonable request.